Trap-induced resonances in controlled collisions of cesium atoms 
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We consider the feasibility of observing a trap-induced resonance [Stock et al, Phys. Rev. Lett. 
91, 183201 (2003)] for the case of two 133 Cs atoms, trapped in separated wells of a polarization- 
gradient optical lattice, and interacting through a multichannel scattering process. Due to the 
anomalously large scattering length of cesium dimers, a strong coupling can occur between vibra- 
tional states of the trap and a weakly bound molecular state that is made resonant by the ac-Stark 
shift of the lattice. We calculate the energy spectrum of the two-atom system as a function of the 
distance between two potential wells by connecting the solutions of the Schrodinger equation for 
the short-range molecular potential to that of the long-range trap in a self-consistent manner. The 
short-range potential is treated through a multichannel pseudopotential, parametrized by the K 
matrix, calculated numerically for atoms in free space in a close-coupling approximation. This cap- 
tures both the bound molecular spectrum as well as the energy-dependent scattering for all partial 
waves. We establish realistic operating conditions under which the trap-induced resonance could be 
observed and show that this strong and coherent interaction could be used as a basis for high-fidelity 
two-qubit quantum logic operations. 

PACS numbers: 34.50.-s, 34.20.Cf, 34.90.+q, 03.67.Lx 
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I. INTRODUCTION 



Ultracold atoms in tightly confining traps are of in- 
terest in a variety of physics disciplines, with important 
applications to the study of quantum degenerate gases 
and the implementation of quantum information proces- 
sors p|. At the heart of these systems are two-body 
atomic interactions whose precise controllability via mag- 
netic Feshbach resonances Q has opened the door to the 
exploration of new phenomena such as the dynamic col- 
lapse of a Bose-Einstein Condensate (BEC) [3j, the pro- 
duction of ultracold molecules Sljj , and the observation 
of the BEC/BCS crossover Strong confinement 

and reduced dimensionality can have a large effect on the 
tunability of the atomic coupling constant 0, 0, 0, 0] . 
For example, quasi one-dimensional (ID) Bose systems, 
where the atomic gas is strongly confined in two direc- 
tions, lead to formation of the Tonks-Girardeu gas whose 
elementary excitations obey Fermi-Dirac statistics |9( . In 
this system, as observed in an optical lattice system of 
tightly confined ID tubes ^3|> a resonance in the effective 
ID coupling strength occurs, leading to strong repulsive 
interactions so that exchange of atoms is suppressed . 
Recently, the origin of this resonance has been attributed 
to a Feshbach-like, confinement-induced resonance that 
is associated with a bound state of a closed scattering 
channel in the transverse mode |14|. This resonance is 
of particular importance since it allows for the control of 
the effective ID interaction strength and the exploration 
of different regimes of the BEC vs the Tonks gas by sim- 
ply varying the confinement of atoms in the transverse 
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direction. 

Concurrently with the exploration of these 
confinement-induced resonances, we have predicted 
a "trap-induced resonance" (TIR) for two atoms 
in separated, tightly confining potential wells, such 
as the sites of a polarization gradient lin-0-lin optical 
lattice |16| . This phenomenon arises when a weakly 
bound molecular state of the dimer is shifted by the 
potential energy (AC Stark shift) of the separated 
wells into resonance with a vibrational eigenstate of the 
trap. In contrast, optical Feshbach resonances |17|. also 
induced by off-resonant laser couplings, arise due to 
strong mixing of excited-state bound molecules with the 
ground unbound states and are generally accompanied 
by inelastic processes. Here, the strong gradients in the 
ac-Stark effect produced by confining optical potentials 
is sufficient to make ground-state molecules resonant 
with atomic trap eigenstates, but with tunability differ- 
ent from the magnetic Feshbach resonances arising from 
the Zeeman shift. The TIR thus provides an alternative 
handle for coherently controlling two-atom ground-state 
interactions. 

In this paper, we evaluate the feasibility of observ- 
ing these resonances in a particularly promising species, 
133 Cs, trapped in optical lattices. The enormous scat- 
tering length of 133 Cs, arising from a very weakly bound 
molecular state near dissociation, leads to the possibility 
of creating trap-induced resonances under realistic ex- 
perimental circumstances. For a complete description of 
the atomic interactions, we develop a multichannel scat- 
tering model based on the generalized pseudopotential 
presented in [l^. Using this model, we calculate the 
two-atom energy spectrum under reasonable operating 
conditions. The remainder of this article is organized as 
follows. In Sec. II we review the basic physics of the 
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TIR and its generalization to a multichannel scattering 
interaction. We apply this formalism to the specific case 
of 133 Cs in Sees. Ill and IV, and conclude in Sec. V. 



II. BACKGROUND 

A. Trap-induced resonances 

We review here the salient features of the TIR 15] . 
Consider two interacting atoms in spatially separated 
harmonic traps, described by the center of mass and rel- 
ative coordinate Hamiltonians, 
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where M = 2m and fi = m/2 are total and reduced 
masses, respectively. Though the two atoms may be 
identical, we assume they are effectively distinguished by 
their internal states. The traps for the two atoms are 
separated in the z-direction by Az. Since the character- 
istic length scale of the van der Waals interaction, /3g, is 
typically much smaller than the size of the trap ground 
state 0] , the interatomic interaction Vi nt (r) can be 
replaced by a contact "pseudopotential" for s-wave scat- 
tering, 



- / x 2nh 2 , . „, N d 

V cS {v,E K ) = a cS {E K )S(r) — r. 
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Note, that we parametrize the pseudopotential with an 
energy-dependent scattering length a e g(E K ) 0, Il9| . 
This is particularly important for strong interactions out- 
side the Wigner threshold regime. In addition, such a 
pseudopotential not only describes the complete energy- 
dependent scattering behavior, but also the bound state 
spectrum if the scattering length is analytically continued 
to negative energies and the Hamiltonian's eigenvalues 
are found self-consistently. 

The energy spectrum associated with Eq. Q is easily 
understood in two limits. When Az — > oo, interactions 
are negligible and the spectrum is a harmonic ladder in 
three dimensions. When Az — > 0, we map to the problem 
of two atoms in a harmonic trap, interacting via a 5- 
function potential in three dimensions, parametrized by 
the scattering length a. The spectrum of this problem 
was solved by Busch et al. l20l| and generalized to higher 
partial wave interactions in Il8t l2lll22( . revealing positive 
(negative) energy shifts of the harmonic ladder for a > 
(<)0 that saturate at Hlu for \a\ — > oo. In addition, for 
a > 0, the spectrum includes a perturbed bound-state of 
the J-function potential. This state represents the near- 
dissociation molecular bound state associated with the 
positive scattering length of an attractive potential. 

For intermediate separations Az, we treat the term 
(jjLijj 2 Az)r cos 9 as a perturbation to the Busch solution 



and rediagonalize the Hamiltonian. This is done for a 
range of scattering lengths to determine E(a). We then 
solve independently the problem of the scattering length 
as a function of kinetic energy, a{Ex) for a given real- 
istic interaction potential between the particles in free 
space. We numerically search for the self-consistent solu- 
tion, setting Ek = E(a) — Vaz, where Vaz = /iw 2 Az 2 /2 
is the potential energy the two particles must overcome 
to collide at zero range. This procedure allows us to sepa- 
rate the long-range behavior of the wave function from its 
short-range behavior and then match them in one self- 
consistent solution. This is particularly important for 
solving realistic atom interaction where molecular poten- 
tials have the range of angstroms whereas the trap scale 
is on the order of many nanometers. 
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FIG. 1: Schematic of trap-induced resonance and eigenspec- 
trum as a function of separation. Fig. (a) shows the sum 
of the harmonic trapping potential and chemical binding po- 
tential (gray line) in the relative coordinate z for zero trap 
separation. Fig. (b). The molecular bound state at Eb and 
trap eigenstate at Strap can become resonant at a critical sep- 
aration Az ros . Fig. (c). Eigenspectrum for two interacting 
atoms in separated traps for a constant scattering length of 
a = lro (ro = y H/Jjujj). Avoided crossings occur in the eigen- 
spectrum as a function of separation Az for all vibrational 
state with zero transverse excitation. 

To verify the use of the energy-dependent pseudopo- 
tential approximation, we have employed a spherically 
symmetric s tep potential as a toy model of the molecu- 
lar binding |15|. For this problem, we can find an an- 
alytic solution for the scattering phase-shift from which 
we deduce the energy-dependent scattering length. The 
self-consistent solution is compared with the full numeri- 
cal solution to the toy problem for a potential with large 
positive scattering lengths at zero energy. For atoms in 
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very tight traps, as the separation is increased, the trap- 
ping potential can be strong enough to raise a bound 
state of the molecular interaction to positive energies [see 
Figs. ma)-njb)]. Avoided crossings occur in the eigen- 
spectrum for certain separations where the energy of this 
molecular bound state becomes resonant with the eigen- 
states of the trap potential [see Fig. dc)]. This effect 
generally occurs if the interaction potential for the free 
atoms possesses a molecular bound state very close to 
dissociation, i.e., the scattering length in free space is 
positive and on the order of or larger than the size of 
the trap ground state. Here, the self-consistent solution, 
which employs an energy-dependent scattering length, is 
crucial; a constant scattering length approximation does 
not capture the location of the bound state and the re- 
sulting energy gap accurately. Also note that in this case, 
analytic continuation of the scattering length to negative 
energies is essential. The TIR occurs through a tunneling 
barrier as the trap distorts the atomic potential at dis- 
tances large compared to the range of the molecular bind- 
ing. As such we probe the "negative energy" spectrum 
of the interaction (with respect to collisions of free parti- 
cles) and the scattering length used for the self-consistent 
solution is taken at the, usually negative, kinetic energy 
in the system. 



B. The multichannel pseudopotential 
approximation 

The model presented thus far involves only single chan- 
nel elastic scattering. A complete analysis of the TIR for 
a given atomic species requires a generalization of the 
pseudopotential to multiple scattering channels, includ- 
ing the full spin-dependent nature of the cold-collision 
process via the hyperfine and exchange interactions, and 
higher-order partial-wave scattering. We thus require 
a multichannel matrix generalization of the scattering 
length that parametrizes the pseudopotential operator, 
thus allowing the inclusion of both elastic and inelastic 
processes. We accomplish this by solving the Schrodinger 
equation for two interacting particles in free space in a 
standard close-coupling expansion of the total wave func- 
tion with good quantum numbers (J, M, IT), the total an- 
gular momentum, its projection on the internuclear axis, 
and parity respectively [23| . 

m, n,E) = J2 (4 M , n ; r *,r) Fl ''f r) ■ (3) 

a 

Here, the indices 7 and 7' refer to the different "in" and 
"out" channels of stationary scattering theory Q ■ We 
choose the channel-state basis to the product basis of 
atomic hyperfine sublevels for the internal degree of free- 
dom, |/i,to/i) I/2, Ti/2}, plus the partial wave quantum 
numbers for the relative motion, \l,m). The functions, 
Fji_ 7 (r), are reduced radial wave functions that satisfy 



the close-coupled equations, 

= £ % (*W'M - K :s r,r) iV.iM . (4) 

7" 

with wave number fc 7 = 2/i(E — E^)/h 2 for the asymp- 
totic channel state energy E^° . In the asymptotic limit 

where Vy, 7 — > 0, the equations decouple and the solu- 
tions are linear combinations of the (reduced) spherical 
Bessel functions defined as J 7 (r) = yHt^rji(k 7 r) and 
Neumann functions Ny(r) = y/k^rni{k^r) , 

F 7 T 7 (r) -f Jy (r)S Yn + N Y (r)K Yn , (5) 

expressed in terms of the K matrix |24j . Typically, 
the scattering length is defined by the scattering phase- 
shift via the S matrix. As we will see, however, in the 
pseudopotential approximation described below, it is ad- 
vantageous and more appropriate to define a scattering 
length matrix through the K matrix. As derived in de- 
tail below, basing the definition of the pseudopotential 
on the K matrix allows us to correctly match the scat- 
tering boundary conditions for the asymptotic wave func- 
tions in Eq. |SJ). In addition, note that our definition of 
the reduced Bessel functions differs in phase convention 
from 23]; for closed channels these functions are purely 
imaginary. 

The objective of the pseudopotential formalism is to 
replace the exact atomic interaction potential at all en- 
ergies, for all partial waves, and for all channels by a 
zero-range pseudopotential. In 0, we showed how such 
a pseudopotential could be constructed for an arbitrary 
partial wave using the limit of a 8 shell of radius s, in the 
limit s — > 0. We generalize here to a multichannel version 
of this pseudopotential, including anisotropic couplings. 
We define an energy-dependent, Z-wave scattering length 
matrix a 7 ^ 7 ', +1 (fc 7 , fey) in terms of the multichannel K 
matrix, 

k^~ I k i ^ ^ (ky,kyi) = K^i^ . (6) 

A pseudopotential with an energy- dependent scattering 
length matrix has the same properties as the single- 
channel pseudopotential, capturing the complete scat- 
tering and bound-state properties of the participating 
channels [TH Hsj . It will further allow us to continue 
the scattering-length matrix element for the channel of 
interest to negative energies, which are crucial for the 
case of atoms in separated traps (see Sec. III). 

Using the definition in Eq. H3 we take as our ansatz for 
the multichannel pseudopotential 

V(r)=J2^m)%,j'(r){l',m'\ (7) 

7,7' 

with 

%,Ar)=cuva% +l5 ^Pi>, (8) 
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where the constants ciji are to be determined. The reg- 
ularization operator, 



(21' + 1)!!(2Z" + 1)!!. In summary, the pseudopotential 
(defined as acting on the reduced wave functions) is 



J+i 



Pi = 



Q21+1 



(21 + 1)! 9r 2Z + J 



(9) V ln ,(r) 



acts near the origin as the identity on regular radial wave 
functions J 7 (r) and the null operator on the i rreg ular 
function N 1 (r). Please note that in contrast to [lg> we 
define here the pseudopotential to act on the reduced ra- 
dial wave functions. The solution to the radial equation 
for the pseudopotential has the form 



- J 7 '(r)5 7 ' i7 + N-y'(r)K^' n , r > s. 

Jy(r)A 7 /, 7 , r < s. (10) 



F>, 

7 ,7 



F< = 

7 i7 



The matrix A 7 ' i7 follows from continuity of the wave 
function across the shell, 

(2Z'-1)!!(2Z' + 1)!! 



(fc r /«)«'+ 1 



Kji 7 . 



(11) 



Inserting these solutions and the ansatz in Eq. JJJ into 
the close-coupled equation Q and leaving the limit s — » 
for later, we need to evaluate the action of the pseudopo- 
tential on the wave function, 



Vfn"(T)F> (r) 



j'+i"+i 6 ( r - s ) 



rl'+l" 



Pi"F>Jr) 



' j'+i"+i S(r - s) 



fc r+i/y+i 



(21" + 1) 



|| ^7":7 



GL2) 



(2Z + 1)!!(2Z' + 1)!! S(r s) d 2l '+* 



(21' + 1)! 



s l + l Q r 2l' + 1 



(15) 

We apply this solution to the case 133 Cs in the next sec- 
tion. 



III. CALCULATION OF THE SCATTERING 
LENGTH MATRIX OF CESIUM FOR POSITIVE 
AND NEGATIVE ENERGIES 

The cold collision properties of 133 Cs are anomalous 
due to the existence of a bound state of the 3 S+ molecu- 
lar potential very close to dissociation [25| and the large 
second-order spin-orbit coupling (2(|. Whereas the for- 
mer property leads to extremely large triplet scattering 
lengths and possible trap-induced resonances, the lat- 
ter implies possible strong inelastic collision channels. 
As such, a realistic modeling of the scattering is nec- 
essary to study the feasibility of observing the TIR. 
We will accomplish this using the numerical solution 
to the close-coupling equation |0J developed at the Na- 
tional Institute for Standards and Technology (23 based 
on the best known cesium dimer interaction potentials, 
recently refined through spectroscopic measurements of 
Feshbach resonances |2y. The multichannel scattering 
solutions will allow us to determine the energy-dependent 
scattering-length matrix as input to the pseudopotential 
in Eq. ©. 



133 



Cs 



having applied the regularization operator near the ori- 
gin. The remaining boundary condition is found by inte- 
grating the close-coupled equations around the shell, 
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FIG. 2: Hyperfine levels of the 6S1/2 ground state of 133 Cs. 
The magnetic sublevels are labeled with letters from a to p. 



Using Eqs. lfTU )l . {IT )) , and 1^2 ) . we find 
(2Z' + 1)!! 



K, 



7 'V', +1/ V+i 

7 



E 

7" 



Cl'V 



j'+r+i i/'+V2 

V,7' S" S r 

s''+'"+ 2 (2Z" + 1)!! 7 "' 7 



(14) 



the 



Since the K matrix is symmetric K~> 7 
close-coupled equation set is fulfilled and using Eq. © 
we find the strength of the pseudopotential is c;/ ; 



Figure [5] shows the atomic hyperfine states and mag- 
netic sublevels in the 133 Cs atomic ground state labeled 
with letters from a to p. The scattering channel of inter- 
est is the two-atom combination in the stretched states, 
\ap), with a total angular momentum projection quan- 
tum number of m to tai — m /i + m f2 = 7. This channel 
is of particular interest since, in the absence of any spin- 
motion coupling, conservation of angular momentum im- 
plies conservation of these quantum numbers. For this 
reason, the equivalent of states a and p for 87 Rb were 
the ones originally considered for entangling atoms via 
ultracold ground-state collisions 28]. For 133 Cs atoms, 
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Channel 


Energy (10 * a.u.) 


atom 1 


m f i 


atom 2 


m j! 
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TABLE I: Participating channel information. The channel 
of interest is \ap) (channel 9), which defines the zero of en- 
ergy (in atomic units a.u.). The first 30 channels are open, of 
which the first eight channels are the ones with both atoms 
in the lower hyperfine manifold. Channels 31 to 49 are closed 
channels with both atoms in the / = 4 manifold. The channel 
states are labeled following the conventions of Fig. [5] (columns 
3 and 5). The table includes the magnetic quantum numbers 
m/ for each atom (columns 4 and 6), the total magnetic quan- 
tum number m tota i (column 7) and the partial wave quantum 
number I (last column). 



this desirable property is not automatically satisfied due 
to the higher-order spin-orbit interaction [2Jj, |26j , which 
couples a multitude of open and closed scattering chan- 
nels. In our calculation, 49 channels were included in the 
close-coupling set (see Tabic Pi. 

The set of close-coupled equations for the reduced ra- 
dial wave functions of the participating channels is inte- 
grated from an initial radius r\ to a final radius r f using 
a standard renormalized Numerov method [2^, E3| for 
both open and closed channels. There, the logarithmic 
derivative matrix m 7 ' i7 (r^) is calculated for the reduced 
wave functions Fy n (r) according to 



(r f ) 



dr-Ff' ,~f{ r )\r=r f 

F Y,-y(r f ) 



(16) 



Matching the logarithmic derivative to the asymptotic 
wave function, we can determine the i-T-scattering matrix 
from m 7 < i7 (ry) according to 

K 1>n = [N^ir^my^r^-N^rf)]- 1 

[J Y (r f )m Yt7 (r f ) - 4,(77)] . (17) 

Here, Jy (rf) and N' ,{rf) are the derivatives of the re- 
duced spherical Bessel functions at r — rf. The K matrix 
then defines the scattering length matrix that is appro- 
priate for use in the generalized multichannel pseudopo- 
tential as discussed in the previous section. Note that 
in order to get a convergent value of the K matrix, the 
final integration radius, rf, can be larger than the typical 
size of the ground state in an optical lattice site. This 
might raise doubts about the validity of the pseudopoten- 
tial approximation in an optical lattice system, since the 
effective scale of the interaction seems to be much larger 
than the trap size. However, as shown in Refs. 0,0], it 
is the characteristic [3$ length scale that affects the valid- 
ity of the pseudopotential approximation. As long as the 
trapping potential does not distort the chemical binding 
potential, it is valid to calculate the energy-dependent 
scattering length in free space and use this as the input to 
the self-consistent solution in the trap. Since (3q — lOOao 
in 133 Cs, and therefore much smaller than typical trap 
sizes in tight optical lattices, we can safely approximate 
the interaction by a pseudopotential. 

As discussed in Sec. lll Al negative energy scattering be- 
comes important for two atoms in nonovcrlapping traps 
due to the tunneling barrier created by the trapping po- 
tential. In this case, the scattering channels of interest 
are closed and we must extend the K matrix to include 
these channels. This can be achieved by simply contin- 
uing the reduced spherical Bessel functions </ 7 (f*/) and 
N 1 {rf) to negative energies by using a purely imaginary 
wave vector, k — in, and allowing a complex argument 
of the Bessel functions (note, in our convention for the 
reduced Bessel function, these functions are then pure 
imaginary). The complete K matrix, which is Hermitian 
for open channels and anti-Hermitian for closed channels, 
is calculated using Eq. \YI\. Figures Et a )"Etf ) show the 
diagonal scattering length matrix element for the \ap) 
channel as a function of the stopping point rf for var- 
ious negative energies. At negative energies, we expect 
the codes to converge in the calculation of a\ ap ) only for a 
small range of r/ . The close-coupled equations need to be 
integrated far enough into the asymptotic regime where 
the interaction potential vanishes. However, for very 
large r, we expect the wave function to diverge exponen- 
tially and the numerical Numerov codes to become un- 
stable. We expect the region of convergence to decrease 
for larger negative energies, since the wave function di- 
verges as exp(+Kr) (see Figs.|3J). As seen in Figs.|3 we 
can reliably calculate the \ap) scattering length matrix 
element for negative energies as low as —10 fiK. The 
calculation errors in the scattering length of Fig. 0] can 
be estimated from Figs. |3]to be around 5% at — 10 /j,K 
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FIG. 3: Calculations of the diagonal scattering length matrix 
element for the \ap) channel as a function of the stopping 
point 77 for different negative energies. 



whereas for smaller negative energies of around — 1 fiK 
the error can be estimated to be smaller than 0.5%. 

A plot of the diagonal scattering-length element as a 
function of energy for the \ap) channel is shown in Fig.0] 
The curve continues smoothly across zero energy and 
also takes a finite value at the bound-state energy as 
expected. At a positive energy, E ks 4 fiK, there is a 
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FIG. 4: Calculations of the diagonal scattering length matrix 
element, derived from the K matrix, for the \ap) channel as 
a function of energy. 

resonance in the scattering length that is associated with 



inelastic processes. Examining the scattering-matrix off- 
diagonal elements, the \ap) \l — 0} channel primarily cou- 
ples to higher partial wave channels in the same hyper- 
fine manifold, channels 10-30 in Tabled and in particular 
to the \ap) \l — 2) and \ao) \l — 2} channels. Exothermic 
couplings that result in both atoms in the lower hyper- 
fine manifold, e.g., couplings to the \aa) \l = 2} channel, 
are one to two orders of magnitude smaller. To better 
characterize the off-diagonal couplings to all channels, 
we will adopt a complex scattering length approach in 
the following. While the K matrix describes the asymp- 
totic boundary conditions and is therefore essential in 
the pseudopotential description, the scattering S matrix 
is a measure of the coupling between different channels. 
These inelastic processes and the resulting loss from the 
channel of interest are therefore more appropriately de- 
scribed by a complex scattering length defined through 
the S matrix, which is related to the K matrix by, 



(18) 



2i\ 



S = [1 + iK] [1 -iK] _1 . 
By unitarity the diagonal S matrix elements are 
S 7n = 1 — Sy, 7 = exp [2i (A + i/j,)] = \S in \ e 

7' 

where we have defined a complex phase shift S = A + i/i. 
This can be used to define a complex energy-dependent 
scattering length [3l| , 



i(3 1 = — 



_ tan (A + ifi) 



kj l+1 



(20) 



The imaginary part of the scattering length is a measure 
of the "loss" from the scattering channel 7 to all other 
open channels. Note that the imaginary part of the scat- 
tering length not only measures the inelastic collisions 
but also includes the coherent couplings to channels in 
the same hyperfine manifold. Figure El shows the real 
and imaginary part of the complex scattering length for 
positive energies. For scattering energies below 4/ii4T, 
the imaginary part is several orders smaller than the real 
part of the scattering length. In this regime, couplings to 
other channels that are due to spin-spin or spin-orbit cou- 
pling can be neglected and the total magnetic quantum 
number mtotal is approximately conserved. This implies 
that, in this energy regime, interactions of atoms in \ap) 
can be treated to good approximation as single channel 
scattering. 

We have also examined the magnetic field dependence 
of the scattering length in the linear Zeeman regime. The 
background scattering length remains large and constant 
over a wide range of magnetic fields and shows three nar- 
row resonances over the calculated range at 11.6, 25, 
and 28 G in Fig. An identification of the different 
resonances, i.e., the labeling of the resonances by the 
corresponding bound-states, is left for future investiga- 
tion. The resonances are very narrow (less than 1 G) 
and should not lead to any strong magnetic field sensi- 
tivity in the atomic interactions. 
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FIG. 5: S matrix based calculations of the real (a) and imag- 
inary part (b) of the scattering length for the \ap) channel as 
a function of energy. 
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FIG. 6: Calculations of the diagonal scattering length ma- 
trix element for the \ap) channel as a function of magnetic 
field. Three narrow scattering resonances can be observed for 
applied magnetic fields in the range from to 50 G. 



IV. TRAP-INDUCED SHAPE RESONANCES 
OF CESIUM IN OPTICAL LATTICES 

We consider a three dimensional (3D) optical lat- 
tice, created by three pairs of linearly polarized counter- 
propagating laser beams. Along one direction we al- 
low for p olarization gradients in a lin-angle-lin configura- 
tion [16j which shifts cr+ and a— standing waves and car- 
ries atoms in appropriate sublevels together from neigh- 
boring sites and probes the TIR. We approximate the op- 
tical lattice wells by harmonic oscillator potentials with 
frequency u> with 



Huj = 2^V PP E R . 



(21) 



V P p is the peak-to-peak depth of the optical lattice po- 
tential and E R is the recoil energy of a 133 Cs atom. 
The Lamb-Dicke parameter, which measures the local- 
ization of atoms in the lattices, n = J Eji/huj = fc/fo, 
where f = y/h/2moj is the full width at half maximum 
(FWHM) of the harmonic-oscillator ground state of a sin- 
gle atom (mass m). As we deal with relative motion of 



two atoms, we choose to scale i n term s of characteristic 
harmonic oscillator units r$ = yjh/ fiu> = ro = 2fo- 

For the observation of a TIR for the |ap)-channel, we 
require ro to be on the order of, or much smaller than, the 
scattering length which takes the value d| ap ) « 2500ao ~ 

132 nm in the negative energy regime of interest. In 
addition, in order to induce a TIR, the trap potential 
at the origin, Vaz = k^- z * \ must be strong enough to 
raise the molecular bound-state at Eb to positive energies 
so that Eb + Vaz > In an optical lattice Vaz is 
limited by V pp , resulting in the constraint ^h\o < Va z < 
V P p. An optical lattice with a Lamb-Dicke parameter of 
7} = 0.25, V pp = 64E R , hj = 16E R , and f w 34 nm 
fulfills all these requirements and still allows a reliable 
harmonic approximation to the optical lattice potential. 
Also, for 133 Cs the (3q parameter is on the order of lOOao 
or 5.3 nm, i.e., much smaller than fo, therefore allowing 
a reliable approximation of the interaction via a (5-shell 
pseudopotential. 

A complete description of interacting atoms in sepa- 
rated traps should include a calculation of the energy 
spectra as a function of trap separation for the scattering 
channel of interest, as well as for all channels with off- 
diagonal couplings to this channel. It furthermore should 
include the state-dependent nature of the optical lattice 
potentials. While this is of great importance for most 
scattering channels, for the \ap) scattering channel, as 
shown in the previous section, the off-diagonal couplings 
can be safely neglected below 4fiK. The calculation of 
the energy spectrum for the \ap) channel using only a 
single channel in Eq. I|15|) with 1,1' = is sufficient as 
long as the appropriate scattering length matrix element 
calculated from the K matrix is used to characterize the 
interaction. 

Figure Ufa) shows the calculated energy spectrum for 

133 Cs atoms in separated harmonic isotropic traps as a 
function of trap separation Az, assuming an optical lat- 
tice with a Lamb-Dicke parameter of r\ = 0.25. The en- 
ergy eigenvalues are calculated self-consistently at each 
separation as described in Sect. Ill Al using the energy- 
dependent scattering length for 133 Cs, which is shown 
in Fig. 0(b) in harmonic oscillator units. This plot of 
the scattering length nicely illustrates the necessity of 
the application of an energy- dependent pseudopotential 
because of the strong energy dependence of the scatter- 
ing length in the energy range of interest. A constant 
scattering-length approximation would misestimate the 
location of the bound state and therefore the size of the 
avoided crossing in the energy spectrum. Figure Eta) 
shows the main signature of the trap-induced resonance 
in the form of the avoided crossing, which in the case of 
133 Cs results in an almost maximal [compare to Fig. 3 
(inset) in Ref. 0|] energy gap of about Huj/2. 

Controlled ultracold collisions, using two-atom Ram- 
sey interferometry in polarization gradient o ptic al lat- 
tices, have been demonstrated by Mandel et al. |32|. Such 
experiments can be used to probe the TIR as shown in 
Fig. |SJ In a linearly polarized lattice, atoms in their 
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FIG. 7: (a) Self-consistent energy eigenvalues as a function 
of trap separation Az calculated using the K matrix based 
scattering length for the |ap)-channel in 133 Cs, which is shown 
in (b) in harmonic oscillator units. 



vibrational ground state are placed in an equal superpo- 
sition of states | a) and \p). These two states are trans- 
ported in opposite directions as the polarizations are ro- 
tated by an angle 9 into a lin-angle-lin configuration. As 
6 approaches 180° (lin-||-lin configuration), neighboring 
atoms in the \ap) channel collide. The collisional inter- 
action results in a relative energy shift for the two atom 
state \ap) = |01) which acquires a phase <f> relative to |00), 
10), and 1 1 1 ) . The interaction entangles these atoms, 
maximally so when <fi — tt. In the experiment of Mandel 
et a/., a filled lattice was prepared in a Mott insulator 
with one atom per site, loaded from a Bose-Einstein con- 
densate [33|. Controlled collisions lead to an "entangled 
chain" of atoms (the ID cluster state [34| in a perfectly 
filled lattice and for tt phase shift). The periodic disap- 
pearance and reappearance of single atom Ramsey fringes 
due to this entanglement is a measurement of the scatter- 
ing phase shift. For the case at hand, assuming adiabatic 
transport, this is the phase shift accumulated along the 
lower energy curve in Fig. |7Ja). At a critical separation 
between the atoms, a large (negative) phase shift will be 
accumulated, providing a signature of the TIR. By con- 
trolling the final separation and hold time of the atoms, 
one can in principle map the full adiabatic energy curve. 



V. SUMMARY 

In this paper we have considered a trap-induced res- 
onance between a molecular bound state and a trap vi- 
brational state for two 133 Cs atoms in overlapping but 
separated traps. To study the interaction between re- 
alistic alkali atoms with hyperfine structure, we gen- 
eralized the (5-shell pseudopotential, introduced in |18j . 
to the case of multichannel scattering parameterized by 
the K matrix. Using a numerical solution to the close- 
coupling equations and a self-consistent solution to the 
Schrodinger equation, we calculated the energy spectrum 
for 133 Cs atoms interacting on the stretched-state chan- 
nel, |/ = 3, to/ = 3) |/ = 4, to/ = 4), as a function of trap 
separation. The large energy gap in the spectrum pro- 
vides a signature by which the trap-induced resonance 
could be experimentally observed for a lattice with Lamb- 
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FIG. 8: Encoding in the 133 Cs hyperfine structure and 
schematic of controlled collisions via trap-induced resonances, 
(a) Hyperfine levels of the 65*1/2 ground state of 133 Cs. The 
logical basis states |0) and |1) are encoded in the stretched 
states in order to avoid inelastic collisions, (b) Following the 
proposal by Jaksch et al. 28] atoms travel on different lattice 
potentials. Unlike in the original proposal, here, the phase 
shift is acquired for separated atoms using the trap-induced 
resonance. 



Dicke parameter ry = 0.25 using Ramsey interferometry 
to probe the collisional phase shift. 

The TIR provides a new method for tuning interactions 
between ultracold atoms, with applications to quantum 
degenerate gases, molecular control, and quantum logic. 
For example, by adiabatically tuning through the TIR to 
zero separation between neighboring traps, atoms in the 
first vibrational state are mapped onto the highest molec- 
ular bound state. In principle, one can induce coherent 
superpositions of atomic and molecular states to control 
new states of matter. The existence of trap-molecule 
couplings opens avenues for new protocols in quantum 
control. 

An important practical consideration, not treated here, 
is the effect of finite temperature. As the TIR is explic- 
itly tied to the vibrational state in the trap, the scattering 
phase shift can be sensitive to motional heating. As seen 
in Fig. ^c), avoided crossings with the molecular bound 
state occur for all vibrational trap states associated with 
zero transverse excitation. For large scattering lengths, 
considered here, the size of the coupling is not very sen- 
sitive to the localization of atoms in the trap and thus 
a common phase shift should be seen for atoms with ei- 
ther zero or one quantum of longitudinal vibration. This 
should make the TIR fairly robust as long as the trans- 
verse excitation is suppressed. An additional consider- 
ation is the effect of spin-dependent trapping. In our 
multichannel scattering calculation, we ignored the spin 
dependent light shift forces that occur in polarization 
gradient optical lattices. Such forces can close certain 
channels and provide new mechanisms for robust encod- 
ings of quantum logic, immune to spin-flips that would 
otherwise occur in free space. These considerations will 
be studied in future research. 
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